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Abstract 

The spin-stiffness of frustrated spin-1/2 Heisenberg models in one and two dimensions is 
computed for the first time by exact diagonalizations on small clusters that implement 
spin-dependent twisted boundary conditions. Finite-size extrapolation to the thermody- 
namic limit yields a value of 0.14 ± 0.01 for the spin-stiffness of the unfrustrated planar 
antiferromagnet. We also present a general discussion of the linear-response theory for 
spin-twists, which ultimately leads to the moment sum-rule. 
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One of the most basic questions in the study of magnetism is that of the existence or 
absence of long-range order in the corresponding magnetic moments. A partial answer to 
this question can be given by the determination of the so-called spin-stiffness of the magnet, 
which measure the rigidity of the spins with respect to a small twist. In particular, systems 
possessing long-range spin-order are stiff, while spin systems that are not stiff accordingly 
show no long-range order in the moments. In the case of spin-1/2 systems, the latter 
stiffness can be directly measured by the generation of a spin-current with a spin-dependent 
magnetic field, as was first shown by Shastry and Sutherland. 1 This method is analogous 
to that used to measure the charge-stiffness of a system, 2 which discriminates between 
metals and insulators. 

In this paper, we apply the former method to the case of the spin-1/2 Heisenberg 
model, Hq = j) ■Aj'^i ' ^j) on both chain and square-lattice geometries. First, we give a 
general discussion of spin-twists for this model based on linear-response theory, which ulti- 
mately leads to the moment sum-rule. 3 We then make practical use of the above method to 
measure the spin-stiffness of near-neighbor Heisenberg chains and square lattices by exact 
diagonalization of the S z = subspace with the Lanczos technique. 4 Employing finite- 
size extrapolations, we find values of the spin-stiffness for the nearest-neighbor Heisenberg 
ferromagnet and antiferromagnet on the square lattice that agree with spin-wave theory 
results to within ten percent. 5 In the particular case of the square-lattice, where reported 
results for this quantity vary widely, 6 we obtain an upper bound for the spin-stiffness of 
p s / J = 0.174, as well as an extrapolated value of p s / J = 0.14±0.01 in the thermodynamic 
limit. Also, for the case of spin-1/2 frustrated antiferromagnets with next (next) nearest- 
neighbor interactions, we generally find that the stiffness coefficient vanishes near the point 
where the analogous classical model losses long-range order in the magnetic moments. 7 We 
now turn to the derivation of the moment sum-rule. 

Linear-Response. Following ref. 1, the rigidity of a Heisenberg model with respect 
to a twist about the spin z-axis is reflected in the ground-state energy of the modified 
Hamiltonian, 



where 6>y represents the twist angle on the bond (i,j). Hence, in the limit of small twists, 
the above Hamiltonian is expressible as H = H Q +Hi, where H represents the unperturbed 
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Heisenberg model with 6>ij = 0, and where the perturbation to this Hamiltonian is given 
by 

#i = 5>4 S) " \%MS?Si + h.c.)]. (2) 

(iJ) 

Above, j i j S ' ) = | Jij(S'. + 5'r — h.c.) is the ^-component of the spin-current operator. Consider 
now the case of identical twists, 9 X , that exist only along nearest- neighbor bonds oriented 
along the x-axis. Then since the spin-rigidity, D s: is related to the difference in the ground- 
state energy by E (9 X ) — E (0) = ND S 9 X for small 9 X , second-order perturbation theory 
gives 



E,, — En 



(3) 



where the spin kinetic energy operator and the spin-current operator along the ^-direction 
are defined by Tj s) = £, \ J iA+ *(SfSr + . + h.c.) and j x s) = £. | J^Sf S{~+* - h.c), 
respectively, and where N denotes the number of sites. 

We can next consider placing a small uniform dynamical twist, 9 x (t) = 9 x e lu3t on 
all of the nearest-neighbor bonds that are oriented along the ^-direction in the modified 
Heisenberg model (1). Application of the Kubo formula then yields a variation in the 
ground state energy per site of N~ 1 AE Q = ^Tl xx 9 x , where 



\{o\A s) W)\ 2 \W\o)\ 



;(*)ln\l2 
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Since the spin-current is given by j x s ^ = dEn/d9 x = NU X S X 9 X , the spin conductivity, 
a s (uj) = Hxx /iu>, is then just 



(4) 



Recr s (o;) =2n 
Im(T s (w) =uj~ 1 
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Integrating Eq. (5a) over all frequencies, and substituting in expression (3) for the spin 
rigidity, we obtain the following moment sum- rule 3 for the spin conductivity: 



duRe <j s (uj) = tv 



/_rp(s)\ 
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N 
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Hence, the fraction of the moment sum-rule occupied by the static twist-response is simply 



where we define p s = 2D S to be the spin-stiffness. Notice that the above result indicates 
that the entire magnetic moment is made up only of excited states in unstiff spin systems 
with p s = 0. Below, the spin-stiffness, as well as the latter static moment-fraction, is 
computed numerically using the Lanczos technique on finite chains and square-lattices for 
both the ferromagnet and for frustrated antiferromagnets. 

Ferromagnet. To check the validity of the method we consider first the nearest- 
neighbor ferromagnetic spin-1/2 chain with periodic boundary conditions. The ground- 
state of Hamiltonian (1) has been obtained by applying the Lanczos technique for N = 
8, 10, 12, 20 sites in the S z = subspace, that permits introduction of twists along the 
z— axis. After finite-size extrapolation to the thermodynamic limit as function of iV -1 , we 
obtain a value for the spin-stiffness of p s j J = 0.248±0.005, which is quite close to the exact 
result of p s / J = s 2 = 1/4 for spin s = 1/2. 1 Similar results were obtained in the case of the 
square lattice. Also, the average spin kinetic energy for the nearest-neighbor ferromagnet 
is simply the total energy in the S z = subspace; i.e, (T^ ) = —Ns 2 J. Hence, Eq. (7) 
indicates that I = 1, which means that the static twist-response saturates the moment 
sum-rule (6) in the case of the nearest-neighbor ferromagnet. We have recovered the latter 
result numerically to within computer accuracy. The saturation effect can also be directly 
understood by the comparison of expression (3) for the spin-rigidity and expression (6) for 
the moment sum-rule, coupled with the observation that the ferromagnetic state is a null 
eigenstate of the spin-current operator. Hence, the spin-response of such a ferromagnet 
is analogous to the charge-response of non-interacting electrons, where the Drude weight 
saturates the / sum-rule. 1 ' 2 

Frustrated Antiferromagnetic Chain. Consider next a periodic spin-1/2 chain with 
both nearest-neighbor (Ji) and next-nearest-neighbor {J 2) antiferromagnetic interactions. 
Again, we have performed exact diagonalizations of Hamiltonian (1) in the S z = subspace 
for N = 8, 10, 12, and 20 sites. The stiffness extracted from these studies are shown 
in Fig.l. After performing a finite-size extrapolation of our results for the unfrustrated 
antiferromagnetic chain (J 2 = 0) as a function of iV -1 , we obtain a value of p 8 /J\ = 
0.270 ± 0.005 for the spin-stiffness in the thermodynamic limit, that is slightly greater 
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than the exact value of 1/4 (see ref. 1). The small discrepancy we obtain with respect to 
the exact result could be due to the fact that the existence of algebraic long-range order 
in the spin-1/2 Heisenberg chain 1 ' 8 exaggerates finite size effects. 9 Also, we see from Fig. 
1 that while the stiffness rises very slightly upon the introduction of frustration, J 2 > 0, it 
drops precipitously to zero around J2/J1 = 0.43. This point is in the vicinity of the well 
studied spin-Peierls dimerization transition, 8 evidenced by the absence of spin-rigidity in 
the chain. Recent estimates that exploit conformal invariance find a value for the latter 
critical frustration of J 2c = 0.24 J 1: 10 which is consistent with the decrease of J 2c with 
increasing lattice size that we observe (see Fig. 1). In fact, finite-size extrapolation of 
these results yields a value of Jicl J\ = 0.33 ± 0.05. We have also computed the fraction, 
Jo, of the moment sum-rule occupied by the static twist-response, which is shown in the 
inset to Fig. 1. In the case of the unfrustrated antiferromagnetic chain (J 2 = 0), this 
value extrapolates to Iq = 0.915 ± 0.005 in the thermodynamic limit, iV _1 — > 0. It is 
intriguing to remark that the latter value is quite close to the analogous fraction of the 
/-sum- rule contributed to by the Drude weight in the t — J model chain with one hole, 11 
which is 0.938. In addition, this fraction increases to a maximum value approaching unity 
at J 2 /Ji = 0.25 of io = 0.986 ± 0.005, just before plummeting to zero. 

Frustrated Antiferromagnet on the Square Lattice. We have also diagonalized Hamil- 
tonian (1) for the case of spin-1/2 on finite square lattices with nearest- neighbor (Ji), 
next-nearest-neighbor ( J 2 ), and next-next-nearest-neighbor (J 3 ) interactions. This model 
has been widely studied because of its close connection with the t — J model on the 
square lattice, 12 and hence because of it's relation to the phenomenon of high-temperature 
superconductivity. 7 ' 13-15 In particular, we have found ground states for N = 8, 16, 18 and 
20 site square lattices with periodic boundary conditions 16 via the Lanczos technique. In 
the case of the nearest-neighbor Heisenberg antiferromagnet we obtain stiffness values of 
Psl Ji = 0.185, 0.177 and 0.174 for systems sizes of iV = 16, 18 and 20, respectively. Notice 
the general downward trend with increasing lattice size. After finite-size extrapolation to 
the thermodynamic limit we find a value of p s j J\ = 0.14 ± 0.01. Both this value and the 
former upper-bound of 0.174 for the spin-stiffness lie below that of p s /Ji = 0.18 obtained 
from second-order spin- wave theory. 5 

Our results for the spin-stiffness of the frustrated J\ — J 2 model (J3 = 0) are shown 
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in Fig. 2, while those of the frustrated J\ — J 3 model (J 2 = 0) are shown in Fig. 3. 
As intuitively expected, we observe that the spin-stiffness generally decreases smoothly as 
frustration increases. In particular, the stiffness vanishes near J2/J1 = 0.5 in the case of 
the Ji — J2 model with 20-sites, whereas it vanishes near J3/J1 = 0.35 in the case of the 
J\ — J3 model. The latter parameter values are close to the points where the corresponding 
classical model looses its long-range Neel order. 7 The fraction, Iq, of the moment sum-rule 
exhausted by the static twist-response for both the J\ — J 2 and the J\ — J 3 models are 
also plotted in the insets of Figs. 2 and 3, respectively. Fig. 2 shows that the J\ — J 2 
model on the 4x4 lattice has a distinct feature near it's classical critical point for values of 
frustration ranging from J2/J1 = 0.55 to J2/J1 = 0.80. Here, the spin-stiffness rises with 
increasing frustration, but then finally vanishes. It has been pointed out in the literature 
that uniform chiral correlations peak near this region, 14 and that the excited states in 
this vicinity are spin singlets. 15 This feature could therefore corresponds to a phase with 
uniform chiral spin order, 13 since the latter state is characterized by spin-0 excitations. A 
spin glass phase, however, is not ruled out. 15 Note that the absence of this feature on the 
other 18 and 20 site lattices that we have studied could be due to their "tilted" nature. 16 
For example, collinear order - which we know must occur in the thermodynamic limit 
for large values J 2 7 - is not possible in such lattices. Clearly, similar diagonalizations of 
Hamiltonian (1) on a 6 x 6 lattice are necessary in order to resolve this issue. 

In summary, we have extended the theory of the determination of spin-rigidity via 
twisted boundary conditions 1 to the general case of quantum Heisenberg models. This 
method has been applied for the first time to the exact diagonalization of frustrated spin- 
1/2 antiferromagnetic chains and square-lattices. Notably, we find an upper bound of 
0.174 for the spin-stiffness of the unfrustrated antiferromagnetic Heisenberg model on the 
square lattice that agrees to within a few percent with spin- wave calculations. 5 However, 
the extrapolated value of 0.14 ±0.01 obtained from finite-size scaling to the thermodynamic 
limit is considerably smaller. 
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Figure Captions 



Fig. 1. Shown is the spin-stiffness, p s , of the frustrated antiferromagnetic chain (in units 
of J\) as a function of next-nearest- neighbor frustration, J2/J1, for various system 
sizes. The inset displays the fraction, I 0: of the sum-rule (6) exhausted by the static 
twist-response. 

Fig. 2. Above, we display the spin-stiffness of the square-lattice Heisenberg antiferromagnet 
(in units of Ji) with only next-nearest- neighbor (diagonal) frustration, J 2 , for various 
system sizes. The inset shows the corresponding fraction (7) of the moment sum-rule 
exhausted by the static twist-response. 

Fig. 3. Similar to Fig. 2, with the exception that only next-next-nearest-neighbor frustration, 
J3, is considered. 
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